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Abstract. We have monitored a Type II outburst of the Be/X-ray binary MXB 0656-072 in a series of pointed RXTE ob- 
servations during October through December 2003. The source spectrum shows a cyclotron resonance scattering feature at 
32.8^ keV, corresponding to a magnetic field strength of 3.67^ JJJ x 10 12 G and is stable through the outburst and over the 



pulsar spin phase. The pulsar, with an average pulse period of 160.4 : 
outburst. From optical data, the source distance is estimated to be 3.9 : 
and a theoretical prediction of the pulsar spin-up during the outburst. 



0.4 s, shows a spin-up of 0.45 s over the duration of the 
0. 1 kpc and this is used to estimate the X-ray luminosity 
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1. Introduction 

MX B 0656-0 7 2 was first classified as a transient source 
bv IClark etaD dl975h when detected at 80mCrab on 1975 
September 20 by SAS-3. Subsequently Ariel V observed the 
source at 50 and 70 mC rab on 1976 March 19 and 27 respec- 
tively ( Kaluzi enskill 19761) . Kaluzienski's comments that MXB 
0656-072 may more closely resemble a long term variable 
source as opposed to a transient possibly led to it bei ng classi- 
fied as a Low Mass X-ray Binary in the catalogue of iLiu et all 

The source, which has been dormant and unobserved 
since 1975/1976, made a reappearance in the X-ray sky in 
a large, extended outburst during 2003 October. During this 
outburst,which lasted over 2 months and reached an X-ray 
luminosity of 200mCrab, it was identified as a pulsating 
X-ray binary with a pulse period of 160.7 s ( M organ et alJ 
2003) an d with an optical counterpart of spectral type 09.7Ve 
llPakulletalJl2003l) . re-categorising it as a High Mass X-ray 
Binary. Time -resolved B and R photometry ( Bart olini et alJ 
2005) of the optical counterpart during this outburst revealed 
no optical periodicity at the X-ray pulse period. Orbital param- 
eters of the system remain undetermined. 

Preliminary ana lysis of the phase averaged spectrum by 
iHeindl et al.1 J2003b showed that the continuum could be fit 
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with a power law with an exponential cutoff. A cyclotron reso- 
nant scatterin g feature with a cent roid energy of 36+1 keV was 
discovered by Heind l et alJ II2003I) . 

The strong magnetic fields found near the polar caps of neu- 
tron stars quantise the ambient electron energies into Landau 
levels. Photons at these resonant electron energies are scattered, 
creating cyclotron resonant scattering features (CRSFs) in the 
observed spectrum. The energy of the fundamental CRSF is 
approximated by: 

where z is the gravitational redshift and B the magnetic field 
strength. Thus measuring the energy of cyclotron features 
in the spectra of accreting X-ray binaries gives us a direct 
measurement of the magnetic field of the neutron star. 



2. Data Reduction 

MXB 0656-072 was observed with RXTE during 2003 
October through 2004 January. Observations coincided with a 
Type II dStellaet alj fl986) X-ray outburst and data were ac- 
cumulated during both the rise and decline of the outburst. In 
total, the PCA livetime for this set of observations amounts 
to 95 ks. Fig. G] shows the MXB 0656-072 lightcurve from 
the All Sky Monitor (ASM) aboard RXTE during the outburst. 
Triangles indicate the dates of pointed RXTE observations. 
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Fig. 1. The ASM lightcurve of MXB 0656-072 during the 2003 outburst with bin sizes of 1 day. Triangles indicate the dates of 
pointed RXTE observations. The arrows indicate the dates at which spectra were compared with each other in order to observe 
any spectral evolution over the outburst. This is further discussed in Sect. 3.2. 



While Standard2f mode PCA data (with 16 s time resolu- 
tion and 128 channel energy resolution) were employed to gen- 
erate phase averaged spectra, the phase resolved spectra and 
timing results were obtained from observation-specific binned- 
mode data with a temporal resolution of 250 ms and an energy 
resolution of 128 channels. 



3. Phase Averaged Analysis 

For the phase averaged spectral analysis, the spectra from 
all observations through the peak of the outburst (from MJD 
52932 to MJD 52964) were added together in an effort to max- 
imise the signal-to-noise ratio. 



Analysis was performed with HEASOFT version 5.3.1 
and spectral fitting with XSPEC version 11.3. It. Phase re- 
solved spectra were created using a modified version of the 
FTOOL fasebin dev eloped at the University of Tubingen 
jKrevkenbohml2004 . 

For spectral fitting, PCA data djahoda et all2 006) in the en- 
ergy range 3-22 keV were used, to which we added systematic 
errors of 0.5%. Below and above this range the systematic er- 
rors required are larger to account for the uncertainties in the 
PCA response matrix, so we ignored data outside the 3-22 ke V 
range. 

F or all spectra we used th e E_8yus_256 JDX1F mode HEXTE 
data (Rothschild et al. 1998), which have a temporal resolution 
of 8 /-is and covers the range 0-250 keV with 256 spectral chan- 
nels. Electronic thresholds limit the HEXTE energy range to 
15-250 keV. In all cases, the signal-to-noise ratio was increased 
by adding together the data from both clusters and generating a 
new response file with weights 1:0.75 to account for the loss of 
signal from one of the detectors in the second HEXTE cluster. 
A binning factor of 3 was applied to data in the 70-200 keV 
range to improve statistical significance of model fitting while 
still maintaining a reasonable energy resolution. Data above 
200 keV were not used. 



3.1. Spectral Model 

In accreting X-ray pulsars the highly ionised accretion stream 
is channelled onto one or both of the magnetic poles of the neu- 
tron star by its strong (B ~ 10 12 G) magnetic field. Depending 
on the accretion rate and the strength of the magnetic field, 
an accretion shock may form. The continuum spectral shape 
is dominated by Comptonisation of photons by the accretion 
stream at the magnetic poles. As a comprehensive theoretical 
spectral model is unavailable, the spectrum is generally mod- 
elled as a exponentially cutoff power law, but see Tanak; I II 19861) 
for the choice of a Fermi-Dirac Cut Off (FDCO) model or 
MihaTal ill 9951) for the NPEX model. In addition to this contin- 
uum we may observe photoelectric absorption by the accretion 
stream or stellar wind, an Fe Ka fluorescence line and one or 
more cyclotron resonance scattering features. 

In our analysis we describe the continuum using a power 
law with photon index F, with a high energy cutoff at E cut . 

I cont {E) = KE r x | l <E - EMa f E I J j (2) 

The energy cutoff is smoothed by including a Gaussian absorp- 
tion line {gsmooth in XSPEC) with an energy dependent width, 
at the continuum cutoff energy. Without a smoothing func- 
tion the energy cutoff forms a discontinuity in the model and 
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Table 1. Spectral parameters (as defined in the text) for the 
phase averaged spectrum. The model consists of power law 
with a smoothed high energy cutoff, photoelectric absorption, 
a CYCLABS cyclotron feature, a Gaussian Fe line and a black- 
body component. The flux in the Fe emission line is repre- 
sented by ape, in photons cm^s" 1 . Uncertainties are at 90% 
confidence intervals. 



Parameter 


Value 




7.28 x 10 21 cm" 2 


r 


89 +0 07 


^eut 


15.5!° 2 keV 


Ef ld 


11.8!° 2 keV 


T 


OMt™* keV 


E c 


32.S + _° 5 4 keV 


w c 


11.8!} ° keV 




0.38^ 2 




6.48+g;g| keV 


0"Fc 


0.363£keV 


fifFe 


0.0195!™°,! 



this results in residuals whic h can easily be mistaken as aris- 
ing fr om a cyclotron feature JKretschmar et aljfl997t Icoburn 
l200lh . Although we attempted to use the FDCO and NPEX 
models to describe the continuum, we found that the best-fit 
continuum was indeed a power law with a smoothed high- 
energy cutoff. 

The spectral shape is insensitive to the H I column density 
and in order to include some estimate of the column density a 
value of 7.28 x 10 21 crrT 2 was esti mated over a cone of radiu s 
1° in the direction of the source (Dickey & Lockm an 119901) . 
The column density, Nh, was fixed at this value for subsequent 
spectral fitting. 

Modifying the continuum is a CRSF modelled with the 
CYCLABS model, which has the analytical form shown be- 
low: 

(W C E/E C ) 2 

CYCLABS (E) = D c ' . (3) 

(E-E c f + Wl 

where E c , D c an d W c are the cyclotron energy, depth an d 
width respectively ( Makishi ma et alJl990tlMihara et al.ll9 90). 
Using other continuum models (FDCO and NPEX) did not al- 
ter the presence of the cyclotron feature. 

We model a strong Fe Ka fluorescence feature with a 
Gaussian emission line and determine the equivalent width to 
be 3.66 keV. In addition we include a blackbody component to 
account for the excess flux at the soft end of the spectrum. The 
blackbody component has a temperature of O^^^keV and 
can be attributed to emission from an accreting polar cap, in a 
manner similar to Cobu rn et all J2001) in the case of X Per. 

The fit to the phase averaged spectrum results in^ 2 ed = 1.1 
with the spectral parameters given in Table ^ The spectrum is 
shown in Fig.|2| 

The cyclotron feature at 32.8^ 4 keV is clearly visible in the 
residuals. A harmonic with width 16! 2 " keV and depth 0. 19^' 2 g 



Table 2. Dates and PCA livetimes of observations used to de- 
termine the spectral evolution over the outburst. 



Date 


MJD 


Livetime (ks) 


2003-10-20 


52933 


11 


2003-10-30 


52943 


8 


2003-11-06 


52950 


15 


2003-1 1-( 16-20) 


5296(0-4) 


7 



may be present at twice the fundamental energy, but detection 
is marginal, as can be interpreted from the uncertainties in the 
line parameters. Assuming that the cyclotron line at 32.8 keV is 
the fundamental (n — 1) and a gravitational redshift of z — 0.3 
for a typical neutron star mass of 1.4M Q and radius of 10 km, 
we can estimate the magnetic field from Eq.l to be 3.67*{Mu x 
10 12 G. 

Note the systematic feature in the residuals at ~10keV. This 
feature is evident in the spectra of a numbe r of ac creting X-ray 
binaries (GS 1843+00, Her X-l) JCoburnlEool. It has been 
observed in spectra accumulated by Ginga llMiharairi995h and 
BeppoSAX (Santangel o et al.ll998h as well, and should be ac- 
counted for in future continuum models of accreting X-ray pul- 
sa rs. A more co mprehensive review of this feature can be found 
in lCoburnl(l200ll) . 

3.2. Spectral Evolution over Outburst 

In order to determine whether there was any significant change 
in the spectral shape or the cyclotron feature during the outburst 
we compared the spectral parameters at various dates through 
the outburst. These dates are indicated by the arrows on Fig.^ 
and further information is shown in Table [2] 

Note that the fourth spectrum is obtained by adding obser- 
vations over a five day range in order to have signal-to-noise 
ratios similar to those of the first 3 datasets. We find that, 
on average, the spectral fits have ;^ 2 ed ~ 1 .2, except for 2003 
November 16-20 which has x^ ed = 1-7. The power law index 
stays between 0.9 and 1 . 1 during this time and there is very lit- 
tle change in E cut and £f id, which have values close to those 
given in Table ^ The cyclotron parameters over this time in- 
terval are shown in Fig. [3] Noting that the error bars represent 
90% uncertainty levels it is clear that the cyclotron feature is 
stable over the outburst. 

4. Phase Resolved Analysis 

If the rotation and magnetic axes of the neutron star are not 
aligned, the observer sees X-ray emission modulated by the 
pulsar spin period. Due to the fact that cyclotron features are 
strongly dependent on viewing angle, phase resolved spec- 
troscopy is a tool well-suited to studying cyclotron features. 

Pulse profiles and phase resolved spectra were generated 
by performing a barycentric correction of the bin times for the 
PCA data and the photon arrival times for the HEXTE data 
for two 15 ks observations at the peak of the outburst. Orbital 
modulation was not accounted for as the orbital ephemeris 
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Fig. 2. The spectrum of MXB 0656-072 from 3-120 keV. The crosses show the data, the smooth curve shows the unfolded 
spectrum and the histogram the model fit. In the second panel the residuals, plotted as Ay 2 , are shown for the case in which no 
CRSF is included in the model. In the lower panel, these residuals are shown again - this time including a CRSF in the model. 
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Fig. 3. Cyclotron line parameters throughout the 2003 October 
outburst. The energy and width are in keV. 



of the system is still unknown. The pulse p eriod was deter- 
mined by running a Lomb-Scarg le analysis JPress & Rvbickil 
1989; Lomb 1976; Scargle 1982) on the barycentre-corrected 
lightcurves. The lightcurves were then folded at the pulse pe- 
riod of 160.4 s to produce the pulse profiles shown in Fig. |4] 

4.1. Pulse profiles 

The pulse profiles are single-peaked and sinusoidal through all 
energy bands and show the pulsed fraction ((F lmx -F m \ n )/F m:ix ) 
increasing steadily from 26% in the 2-5 keV band up to 40% 
in the 30-60 keV band. The pulse shape is consistent with what 
one would expect from viewing just a single magnetic pole ra- 



diating isotropically as a blackbody, including the effec ts of 
gravi tational bending around the neutron star jBeloborodovl 
2002). The almost perfectly sinusoidal pulse shapes may indi- 
cate that there is no absorption along the accretion column, as 
it is believed that strong absorption along the accretion column 
contribut es to the c omplexity o f the pulse profile a t low ener- 
gies JKrevkenbohm et all2002HNagase et alJl983h . The pulse 
profi le shape and a mplitude support the conclusion reached by 
Buli k et alJ (J2003): that the rotation and magnetic axes of ac- 
creting pulsars are almost aligned. 

4.2. Spectral evolution over Pulse 

The data were split into 5 phase bins for analysis, as shown 
in the topmost panel of Fig.|4] This gives us good statistics at 
a temporal resolution sufficient to describe this particular pulse 
shape. Evolution of the spectral parameters with phase is shown 
in Fig.|5] 

For the phase resolved spectra, we used the same model 
as in Sect. 3.1, however the addition of a blackbody compo- 
nent did not significantly improve the fit, so was not included. 
Furthermore, a Gaussian absorption component at ~10keV 
was required to take account of the systematic feature at this 
wavelength, mentioned in Sect. 3.1. 

The cyclotron feature at energy around 32 keV is present 
through all pulse phases. There are, however, variations in the 
cyclotron line width. In general, the line is at its narrowest 
(6*g keV) in the peak rise and reaches its broadest (14^keV) 
during the pulse decline (See Fig.|5]) The cyclotron line depth 
also shows some variation through phase, in a manner in- 
verse to the continuum power law normalisation coefficient. 
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Fig. 4. Pulse profile shown in counts s Profiles are plotted 
twice for clarity. 



However, the cyclotron feature in this source shows little vari- 
ability with phase in comparison with som e other sources 
such as Vela X-l jKrevkenbohmetalJl2002l) and GX 301-2 
jRrevkenbohm et alJl2004l) . which show cyclotron line energy 
variations of a couple keV throughout the phase. 



5. Discussion 

5.1. Source Distance 

As the optical counterpart to MXB 0656-072 has been identi- 
fied, we can use the optical magnitude and colour to estimate 
the distance to the source. 

The optical counterpart has been identified as an 09.7Ve 
spectral type star, with an obse rved B - V colour in the range 
0.86 to 1.02 JPakull et all2003l) and an observed V-magnitude 
in the range 12.05 mag to 12.38 mag. For this estimate we 
shall assume that the variability in the V-magnitude and the 
B - V colour is due to the circumstellar disk contributing the 
the V-magnitude. Hence we will use the faintest available V- 
magnitude value, assuming that this is when the disk is small- 
est, along with the corresponding bluest B - V colour. We have 
estimated errors typical for photometry of a 12th magnitude 
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Fig. 5. From top to bottom the power law index, the power law 
normalisation, the cyclotron energy, cyclotron line width and 
cyclotron line depth are plotted as a function of phase. Error 
bars indicate 90% uncertainty levels. 

star on a 1-2 m telescope and use V = 12.38 + 0.02 mag and 
B-V = 0.86 + 0.03. 

Interstellar reddening is given by: 



E(B -V)-(B-V)-(B- V) 



(4) 



where (B - V)o is the intrinsic colour. T hus for an 09.7Ve star 
we expect (B - V) = -0.28 jWegnerll l994) and an absolute 
magnitude of My — -4.2 mag (Zombeck 1990). 

This puts the magnitude of reddening in the V-band at A — 
3.7 + 0.1 mag. Knowing the apparent magnitude, the absolute 
magnitude and the reddening magnitude allows us to estimate 
the distance using: 



D = 10 



l+(m-M-A)/5 



= 3.9 + 0.1kpc 



(5) 



It is most likely that the optical colour and magnitude used 
in this calculation were observed at time when the contribu- 
tion from the Be star disk was not negligible. Hence the above 
distance estimate will be an underestimate. 

5.2. Spin-up trends 

Accretion of material onto a pulsar is expected to transfer a 
significant amount of angular momentum to the neutron star. 
Thus an accreting X-ray pulsar undergoing an outburst is in a 
state where the mass (and thus momentum) transfer rate is high 
and one expects to observe a slight decrease in the pulse period 
as it is spun-up. 

We have analysed all pointed PC A data in the 2-20 keV 
energy range and determined pulse periods as a function of 
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time. Initially, to deri ve the pulse perio d, we used a Lomb- 
Scargle analysis jPress & Rvbickil ll989l) . Since the resulting 
values h ad a rather high uncertainty, the method of phase con- 
nection JPeeter et alJll98lt iMunoet alJ [2002) was applied in 
a second step: average pulse profiles (with an initial period of 
P = 160.4 s) for data sets of individual RXTE orbits (of typ- 
ical length 1 to 3 ks) were produced and phase connected by 
an integer number of pulse cycles. Using an improved solution 
for the pulse ephemeris, including P, also the larger gaps be- 
tween data sets could be uniquely bridged. The result is shown 
in Fig. [6] The measurements are consistent with an initial pe- 
riod of 160.58 s at MJD 52930 and an average spin-up until 
MJD 52971 of P = -0.0101 + 0.0003 s day" 1 . 

After MJD 52971 the accretion rate drops to a level low 
enough to halt spin-up of the neutron star. This is evidenced by 
a flattening of the slope in the pulse period vs time plot shown 
in Fig. [6] The data after MJD 52971 are consistent with a con- 
stant period. This levelling off of the pulse period to a constant 
value after MJD 52971 indicates that the pulse period changes 
are an effect of accretion torques rather than modulation by the 
orbital period of the neuton star. 

From the average X-ray flux in the 2-10keV band (4.35 x 
10~ 9 erg cirT 2 s -1 ) we can calculate the X-ray luminosity if 
the distance to the source is known. Using our estimate of 
3.9 + 0.1 kpc for the source distance from Sect. 5.1, the X-ray 
luminosity of the source is 6.6 + 0.4 x 10 36 erg s _1 . Then a 
straightforward applic ation of the accretion torque model by 
iGhosh & Lambllfl9 79) (Eq. (15) with standard neutron star pa- 
rameters and our own values of B = 3.67 x 10 12 G and L x = 
6.6x 10 36 erg s -1 ) yields a spin-up rate of P = -0.0034 s day -1 . 
This is about a factor of three smaller than the observed value. 
We do not consider this a problem, as accreting binary pulsars 
very generally show a broad range of behaviour with respect to 
the develop ment o f thei r spin periods with time (for reviews 
see iNagasel Jl989h and Bildsten etal] il997l) '). This includes 
episodes of spin-up and spin-down, depending on the exact 
conditions of accretion at the time. For a transient in outburst 
such as MXB 0656-072 one would expect to observe a strong 
spin-up, since under the greatly increased mass accretion rate, 
as evidenced by the large increase in X-ray luminosity, the sys- 
tem will tend to be far from any equilibrium condition. 



6. Summary 

From analysis of the pointed RXTE observations of MXB 
0656-072 and limited optical data available we can deduce the 
following: 

- The X-ray spectrum continuum is best fit by a power law 
model with an exponential cutoff, where a Gaussian func- 
tion with an energy dependent width is used to smooth dis- 
continuities at the cutoff energy . 

- A cyclotron resonance scattering feature is present at an 
energy of 32.8^1 keV and the energy is independent of the 
continuum model used to describe the spectrum. The CRSF 
is present at this energy in all observations throughout the 
outburst. The CRSF is present throughout the entire neutron 




52940 52960 52980 53000 53020 53040 
JD - 2400000.5 



Fig. 6. Pulse period as a function of time through the outburst, 
with a linear fit to the data. 



star spin phase, but the width and depth of the feature show 
variation through the spin phase. 

- Although there is no convincing evidence for higher har- 
monic cyclotron lines in this data, the possibility of a higher 
harmonic cannot be ruled out. 

- The magnetic field inferred from the cyclotron energy (by 
assuming that the detected line is the fundamental CRSF) 
is 3.67+^ x 10 12 G. 

- On the basis of optical photometric data, the distance to 
MXB 0656-072 is estimated to be 3.9 + 0.1 kpc. 

- Over the time considered, we have measured a strong spin- 
up (about 0.45 s in 30 days), which appears quite reason- 
able under the conditions observed in this high mass X-ray 
binary in outburst. 
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